
> library(foreign)
> x <- read.spss("E:\\EDA FOR ZIPPING\\DEMO EDA DATA\\beck.sav")
> x1 <- data.frame(x)

> 
> #install.packages("moments","e1071")
> 
> library(moments)
> kurtosis(x1$beck,na.rm=T)
[1] 6.599262
> skewness(x1$beck,na.rm=T)
[1] 1.497913
> 
> library(e1071)                    

Attaching package: e1071

The following objects are masked from package:moments:

    kurtosis, moment, skewness

Warning message:
package e1071 was built under R version 3.6.3 
> skewness(x1$beck,na.rm=T)  
[1] 1.496713

# note there are two ways of defining kurtosis with one three less than the other.
# which is excess kurtosis over that expected from a Normal distribution 
# since a Normal distribution has a kurtosis of 3.
# the kurtosis below (3.59) is what we compare to 3 and the above one to 6
# see https://imaging.mrc-cbu.cam.ac.uk/statswiki/FAQ/Simon.

> kurtosis(x1$beck,na.rm=T)
[1] 3.592214
> 
> wilcox.test(beck ~ sex)

        Wilcoxon rank sum test with continuity correction

data:  beck by sex
W = 432582, p-value = 2.108e-07
alternative hypothesis: true location shift is not equal to 0
